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We calculate the Collins fragmentation function in the framework of a spectator model with 
pseudoscalar pion-quark coupling and a Gaussian form factor at the vertex. We determine the model 
I"""- parameters by fitting the unpolarized fragmentation function for pions and kaons. We show that the 

Collins function for the pions in this model is in reasonable agreement with recent parametrizations 
obtained by fits of the available data. In addition, we compute for the first time the Collins function 
for the kaons. 
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i 1 I. INTRODUCTION 

p , 

The Collins fragmentation function [l| measures how the orientation of the quark spin influences the direction of 
emission of hadrons in the fragmentation process and can thus be used as a quark spin analyzer. It contributes to 
several single spin asymmetries (SSA) in hard processes, such as semi inclusive deep inelastic scattering (SIDIS), pp 
collisions and e + e~ annihilation into hadrons. We shall henceforth use the term "Collins asymmetries" to denote any 
asymmetry where the Collins function plays a role. 

The first experimental evidence of a nonzero Collins function for pions came from the measurement of a Collins 
asymmetry in SIDIS on a proton target by the HERMES collaboration [p . The same asymmetry, but on a deuteron 

j. target, was found to be consistent with zero by the COMPASS collaboration Q. At the moment, the most convincing 

(f-) , evidence of a nonzero pion Collins function comes from the measurements of a Collins asymmetry in e + e~ annihila- 
(T) ■ tion [J]. First extractions of the pion Collins function were performed in Ref. [1, @. A recent fit to SIDIS and e + e~ 
annihilation allowed the simultaneous extraction of the Collins fragmentation function and of the transversity parton 
distribution function Q , clearly showing the importance of the Collins function as a tool to investigate the structure 
| of hadrons. The kaon Collins function is at the moment unknown. 

A few model calculations of the Collins function for pions have been presented in the literature [1, [H, E3, EH E2 
and used to make predictions and/or compare to available data EH EH EH ■ However, the above calculations have 
been found to be inadequate to describe the data. 

The aim of the present work is to show that a Collins function in reasonable agreement with the available 
parametrizations can be obtained in a model with pseudoscalar pion-quark coupling and Gaussian form factors at 
the pion-quark vertex. We also present, for the first time, the Collins function for the fragmentation of quarks into 
kaons. This calculation is relevant for the interpretation of recent kaon measurements done at HERMES [lj| as well 
as COMPASS 1171 and for future measurements at BELLE and JLab. 
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FIG. 1: Tree-level diagram for quark to meson fragmentation process. 



II. MODEL CALCULATION OF THE UNPOLARIZED FRAGMENTATION FUNCTION 



In the fragmentation process, the probability to produce hadron h from a transversely polarized quark q, in, e.g., 
the qq rest frame if the fragmentation takes place in e + e~ annihilation, is given by (see, e.g., [Hj]) 

D h/q , (z, Kl) = Df(z, Kl) + Hi\z, Kl) {kX ^l' Sq , (1) 

where Mh the hadron mass, k is the momentum of the quark, s q its spin vector, z is the light-cone momentum fraction 
of the hadron with respect to the fragmenting quark, and the component of the hadron's momentum transverse to 
k. D\ is the unintegrated unpolarized fragmentation function, while H^~ q is the Collins function. Therefore, H^~ q > 
corresponds to a preference of the hadron to move to the left if the quark is moving away from the observer and the 
quark spin is pointing upwards. 

In accordance with factorization, fragmentation functions can be calculated from the correlation function [l9| 



A(z,k T ) = ^- Jdk+A(k,P h ) 



1 ^ [ d£ + d 2 £ T 



2z ^ 

x 



(2) 



with k~ = Pf~ jz. A discussion on the structure of the Wilson lines, U. can be found in Ref. Here, we limit 

ourselves to recalling that in Refs. [13, HH it was shown that the fragmentation correlators are the same in both semi- 
inclusive DIS and e + e~ annihilation, as was also observed earlier in the context of a specific model calculation [2(| 
similar to the one under consideration here. In the rest of the article we shall utilize the Feynman gauge, in which 
transverse gauge links at infinity give no contribution and can be neglected 0, [H, H3 • 

The tree- level diagram describing the fragmentation of a virtual (timclike) quark into a pion/kaon is shown in 
Fig. Q] In the model used here, the final state \h, X) is described by the detected pion/kaon and an on-shell spectator, 
with the quantum numbers of a quark and with mass m s . We take a pscudoscalar pion-quark coupling of the 
form q n -nl T>Ti. where Ti are the generators of the SU(3) flavor group. Our model is similar to the ones used in, e.g., 
Refs. |25l. 12a. I27L l28|. The most important difference from previous calculations that included also the Collins function, 
i.e., those in Refs7[8|, 0, [l^, [ll|, [l2j, is that the mass of the spectator m s is not constrained to be equal to the mass 
of the fragmenting quark. 

The fragmentation correlator at tree level, for the case u — > 7r + , is 

A {0) (k, P ) = 75 - f h + m sh5 |±^1 2nS((k - P h f - ^) (3) 

and, using the (5-function to perform the k + integration, 

A(0)(Z ' M = 32^ (l-z)P fc -(fc»-m^ ' (4) 

where k 2 is related to k\ through the relation 

k 2 = zk\j{\ - z) + m 2 J{\ - z) + M?Jz (5) 
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which follows from the on-mass-shcll condition of the spectator quark of mass m s . We take m to be the same for u 
and d quarks, but different for s quarks. Isospin and charge-conjugation relations imply 

= D d ^ + = Df**~ = Df^*~ , (6) 
D{-* K+ = D\^ K ~ . (8) 





= Df^ + 


= Df~^ w ~ 


= Dp 1 " 




= D{^ + 


= Df~* w ~ 




D d r K+ 


= D d r K+ 
Dr K+ 


= D d r K ~ 
= D s r K ~ 


= D d r K 



For later purposes it is useful to spell out the relations coming from isospin and charge- conjugation relations for the 
unfavored functions 

(9) 
(10) 

(12) 

We assume the above relations hold for all fragmentation functions, in particular for the Collins function. 
The unpolarized fragmentation function is projected from Eq. (j4|) 

D^z, k 2 ) = Tr [A (z, fc T )7 + ]/2 (13) 

leading to the result 

a 9 2 q7I [z 2 k 2 T + (an + m s - mf] 

D ^T)-^ z3 (fc 2 +£2)2 , (14) 

with 

, (1 — z) o n rrvL — m 2 
L 2 = V M 2 + m H . 15 

z 2 z 

In the limit m s = m and setting the form factor to 1, our result for D\ reduces to Eq. (3) of Ref. [HJ (multiplied by 
two because in that paper the results refer to in 7r°). The two nonzero kaon fragmentation functions D\~^ K and 
jjs-^K + are gj ven by tjjg same functional form, but with different masses m, m s , Mh- 
The integrated unpolarized fragmentation function is defined as 

/>OC 

D 1 (z) = nz 2 dk^Di(z,k^). (16) 
Jo 

Here the integration is over the transverse momentum of the produced hadron Kt = —zkx with respect to the quark 
direction, which is why an extra factor of z 2 appears in the above equation. The above integral is divergent. In 
Ref. [lH, a cutoff on kr has been used. On the other hand, in Ref. [13], a Gaussian form factor depending on k\ has 
been introduced at the pion-quark vertex, which effectively cuts off the higher kx region in the integration. Similarly, 
in this work we use a Gaussian form factor of the form 

ggw ^iTL e -^ with A 2 = \ 2 z a (l-zf. (17) 

at the pion-quark vertex. Due to Eq. ([5|), the above form factor effectively cuts off the higher kr region of the 
integration. The form of the vertex is chosen merely on the basis of phcnomcnological motivations, in order to 
reasonably reproduce the unpolarized fragmentation function. With the chosen form factor, the integration in Eq. (|16|) 
can be carried out analytically and yields: 



(1 - z)((m B - m) 2 - Aljty cxp(- ^ ) 



+ [z A A J -2zl(m.-my-M£) \ - H T 0, 



(l-z)A^ 
L 2 „s 2zL 2 



A 2 V ' (1 - z)A 2 
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where the incomplete gamma function is, 



T(0,z) = / — dt. 



(19) 



The parameters of the model are A, a, (3, together with the mass of the spectator m s and the mass of the initial 
quark m. For the latter, we choose a constituent quark mass m = 0.3 GeV for the u and d quarks, and m = 0.5 
GeV for the s quark. To fix the values of the other parameters, we performed a fit to the parametrization of Ref. [29| 
(NLO set) at the lowest possible scale, i.e., Qo = 0.4 GeV 2 . The resulting values for the parameters are 



g qv =4.78, 



A = 3.33 GeV, 



a = 0.5 (fixed), 



/3 = (fixed), 



(20) 



which are common to both pion and kaon fragmentation functions. The only parameters that change according to 
the type of fragmentation function are 



U , — > 7T 

u -> K + 
s -► K + 



m s = 0.792 GeV, 
m s = 1.12 GeV, 
m, = 0.559 GeV, 



m =0.3 GeV (fixed), 
m =0.3 GeV (fixed), 
m =0.5 GeV (fixed). 



(21) 
(22) 
(23) 



Obviously, also the mass of the hadron changes: we take m/, = 0.135 GeV for the pions and irih = 0.494 GeV for 
the kaons. We remark that it is not possible to estimate the errors in the parameters in a meaningful way because 
the fragmentation functions in Ref. 29j have no error bands. It could be in principle possible to use the recent 
parametrizations with error bands [30j, but the lowest scale they reach is 1 GeV 2 , which we consider to be too high 
to compare to our model. 



Fig. [5] show the plots of the unpolarized fragmentation function D\(z) multiplied by z for u 
s — > K + . The parametrization of [29[ (NLO set, Qo = 0.4 GeV 2 ) is also shown for comparison. 



K + , and 



III. MODEL CALCULATION OF THE COLLINS FRAGMENTATION FUNCTION 

We use the following definition of the Collins function [l^ 1 

Fj Ht(z, fc 2 ) = I Tr[A(z, (24) 



M h 1 v 1 11 2 

As is well known [l2||, using the tree- level calculation of the correlator function is not sufficient to produce a non- 
vanishing Collins function, due to the lack of imaginary parts in the scattering amplitude. In order to obtain the 
necessary imaginary part, we take into account gluon loops. In fact, gluon exchange is essential to ensure color gauge 



1 The factor 1/2 is due to a slightly different definition of the correlator in Eq. J2j with respect to Ref. [l2ll 
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FIG. 3: Single gluon-loop corrections to the fragmentation of a quark into a pion contributing to the Collins function in the 
eikonal approximation. "H.c." stands for the hermitian conjugate diagrams which are not shown. 



invariance of the fragmentation functions. Contributions come from the four diagrams in Fig. [3j Diagrams (a) and 
(b) represent the quark self-energy and vertex diagrams, respectively. Diagrams (c) and (d) can be called hard-vertex 
and box diagrams, respectively. For the calculation of the diagrams with the eikonal line, the Feynman rules to be 
used are 



= ig t a K , 



/ 



-l~±it 



(25) 



where g is the QCD coupling. Note that the sign of the ie in the eikonal propagator is different for semi-inclusive 
deep inelastic scattering (+) and e + e~ annihilation (— ), but this does not influence the computation of the Collins 
function [2(| Hl[ . The resulting formula for the fragmentation correlator corresponding to Fig. [3] (d) is 



4a s iQt + in) 



9 q ^l5 (fi - fh + m a ) 2ir S((k - P h f - ml) 



(2tt) 3 k 2 -m 2 

dH i^t a jjjk - f h - f+ m s ) g^ 7 5 iiP - f+ m) i {-ig"-) jit a ) 

(2tt) 4 ((jfe - P h - I) 2 -m 2 s + ie)((k - I) 2 - m 2 + ie)(-l~ ± ie)(l 2 - m 2 + ie) ' 



(26) 



Note that one of the form factors g q7T is inside the loop integral. When using a form factor as in Eq. (JT7J) , it would 
seem reasonable to replace k 2 with (k — I) 2 . However, since the form factor is introduced to the purpose of cutting 
off the high-^T region, we prefer to maintain the form factor depending only on fc 2 , so that it can be pulled out 
of the inte gral and simplify the calculation. This choice is similar to imposing a sharp cutoff on k 2 — as done in 
Rcf. HHlfljj] — and not on (jfe - I) 2 . 

The Collins function is given by (we take always m g = 0) 



Ht(z,k 2 T ) 



2 a s gL 



9" r< 



M h 



1 



HtUz, k 2 T ) + Hh h) {z, k 2 T ) + Ht (d) (z, k 2 T 



(27) 



(2tt) 4 z 2 (1 - z) k 2 

where the subscripts in the r.h.s. refer to the contributions from diagrams [3] (a), (b) and (d) plus their Hcrmitean 
conjugate, respectively. Diagram (c) gives no contribution to the Collins function. The separate contributions read, 
for the fragmentation of u — > 7r + , 

,2^ 



Hu a )(z, k T ) 



k 2 



m 



fig, 



(28) 



Hi( b )(z,k^) = 2m s I 2g , 



Hi( d ){z, k%) = { _/ 34 3 {2zm + m s - m) + I 2g 2zm(k 



m 2 +M 2 (l-2/z) 



+ (m s - m)((2z - l)k 2 - M 2 + m 2 - 2zm{m + m s )j |. 



(29) 
(30) 
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FIG. 4: Half moment of the Collins function for u — > ir + in our model, (a) H^ 1 ^ 2 ^ at the model scale (solid line) and at a 
different scale under the assumption in Eq. (|37[1 (dot-dashed line), compared with the error band from the extraction of Ref. [f|, 

(b) H^ 1 ^ 2 ^ / Di at the model scale (solid line) and at two other scales (dashed and dot-dashed lines) under the assumption in 
Eq. (|38[) . The error band from the extraction of Ref. 0] is shown for comparison. 



The loop integrals I\ g , l2 g and l34 g [12|] 2 are given by 



'is 



l 2g 



7r v /A(m, m g ) 



2k 2 

7T 



hi 



3ig 



2y/X(m a ,M h ) 

F(i - z) 



tt In 



k 2 



k 2 + m 2 s - Ml - y/X(m a ,M h ) 



M 2 + yJ\{m Sl M h ) 



(31) 
(32) 

(33) 



where A (7711,7712) = (k 2 — (mi + 7712) 2 ) (k 2 — (mi — 77i2) 2 )- In the limit rn s = m and setting the form factor to 1, our 
result for H ^ reduces to Eq. (15) of Ref. [l2| (multiplied by two because in that paper the results refer to u — > n°). 
It is important to note that the Collins function should obey the positivity bound [9|, [31j 



Hj-{z,k 2 T ) < Z)i(^,4) 



Integration over k 2 -, gives the simplified expression 



±(1/2) 



(z) 



D 1 (z) 



1 

^2' 



where the half moment of the Collins function is defined as 



H 



±(1/2) 



(z) 



2M h 



Ht(z,k 2 T ). 



(34) 



(35) 



(36) 



In Fig. [4] (a), we have plotted the half moment of the Collins functions vs. z for the case u — > n + . In the same panel, 
we plotted the 1-cr error band of the Collins function extracted in Ref. 0] from BELLE data, collected at a scale 
Q 2 = (10. 52) 2 GeV 2 . In order to achieve a reasonable agreement with the phenomenology, we choose a value of the 
strong coupling constant a s = 0.2. Such a value is particularly small, especially when considering that our model has 
been tuned to fit the function D\ at a scale Qq = 0.4 GeV 2 , where standard NLO calculations give a s rts 0.57 {29l I32I] . 
In any case, the problem of the choice of a s is intimately related with the problem of the evolution of the Collins 
function (see below). 



2 The expression I^g used here corresponds to Eq. (24) in Ref. [l2ll multiplied by k 
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FIG. 5: Half moment of the Collins function for u — > A" + in our model, (a) H^ {1/2) at the model scale of 0.4 GeV 2 , (b) 
jy^-C 1 / 2 ) a t the model scale (solid line) and at two other scales (dashed and dot-dashed lines) under the assumption in 
Eq. J3S). 



In Fig.[4](b), we have plotted the ratio H^ 1 ^ 2 ^ / Di and compared it to the error bands of the extraction in Ref. 0]. 
Also in this case the agreement is good, with the above mentioned choice of a s = 0.2. 

At this point, some comments are in order concerning the evolution of the Collins function (or of its half- moment) 
with the energy scale. Such evolution is presently unknown, except for some work done in Ref. (33| , which is however 
based on questionable assumptions. Some authors (e.g. Refs. @, 0]) assume 



H 



±(1/2) 



Di 



H 



±(1/2) 



Di 



(37) 



Q- 



i.e., that the evolution of H^ 1 ^ 2 ^ is equal to that of D±. This seems unlikely, in view of the fact that the Collins 
function is chiral-odd and thus evolves as a non-singlet. An alternative choice could be to assume 



H 



±(1/2) 



Ql 



= H^ 1/2) 



(38) 



i.e., that H^ 1 ^ 2 ^ docs not evolve with the energy scale. This is an extreme hypothesis, which cannot be true because 
at some point the positivity bound ([33)) would be violated at large z. We demonstrate this in Fig.0](b) where we show 
how the ratio H^ 1 ^ 2 ^ / D\ behaves at at three different energy scales if only D\ is evolved (we use the unpolarized 
fragmentation function of Ref. (29| for this purpose). Clearly, in this case the ratio grows more steeply with z at higher 
energies, due to the decreasing of D\ in the large-z region. While the evolution of the T-odd parton distribution and 
fragmentation functions remain an outstanding issue, these results show that different assumptions on the Collins 
function scale dependence have a significant impact and should be considered with care. 

For the fragmentation u — > K + ', the same analytic formulas are used but with the other set of parameter values. 



IV. ASYMMETRIES IN e+e" ANNIHILATION 



The BELLE collaboration has reported measurements of various asymmetries in e + + e — > 7r^ + ^ + X that can 
isolate the Collins functions Q. In particular, the number of pions in this case has an azimuthal dependence [3~i| 



N hl h 2 (zi,z 2 ) cc^2e q 2 I D 1{q ^ hl) (z 1 )D 1{ g^ h2) {z 2 ) 
q V 



v (39) 



sin 2 9 



where 12 are the azimuthal angles of the two pions relative to their jet axes (or thrust direction) and the 2 jet 
production plane. Normalizing this distribution and extracting the azimuthal asymmetry gives a measure of the 
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FIG. 6: Half moment of the Collins function for s — » K + in our model, (a) tf 1 L(1/2) at the model scale of 0.4 GeV 2 , (b) 
Tjj 1 ' 1 / 2 ) J ]j 1 a t the model scale (solid line) and at two other scales (dashed and dot-dashed lines) under the assumption in 
Eq. J3S). 



product of moments of Collins functions. BELLE noted that there are QCD radiative corrections that compete with 
the leading twist effects. To cancel out those corrections they take the ratio of the asymmetry for unlike sign events 
(it + it~) to the asymmetry for like sign events. This super ratio has the form Q 

*M^) = ^f%(Pu-P L ), (40) 

where 

V ^ ( U Ul '^ \ TIM 1 / 2 ) (, \ i ljMl/2) , \ rr-L(l/2) , n\ 

nj = t r , (41) 

S g C q [ D l (g-Mr+)(*l) D l (q^-){z 2 ) + D l (*! X ) D l (q^+) faj) J 

^ 2 /'h- L ( 1 /2) s rr-L(l/2) / \ i tt-X(1/2) / \ rr-L(l/2) , A 

p _ L ? e g ( g l(g^+)( 2 l) g l( 9 -->.+ )( Z 2) + g 1( ^-)(^l)g 1( ^-)U2)j 

"L — / \ • (42) 

S g e ? ( k - D l(g^7T+)(2l)-Dl(«^r+)(Z2) + A (g-^n-) («l) -^1 (g-^n-) (^) J 

Note that Eq. (|40p is a linear approximation for << 1. For numerical studies, we take the unpolarizcd fragmenta- 
tion functions from Ref. [1^] (NLO set) at the scale of the BELLE measurements, i.e., Q 2 = (10. 52) 2 GeV 2 . We take 
also (sin 2 9) /(l + cos 2 9) « 0.79. 

For the calculation of the asymmetry we have to make some assumptions on the unfavored Collins fragmentation 
functions. In order to have a guiding principle for our assumptions, we consider the Schafer-Teryaev sum rule [35| . 
which states that 



£jf dzH^l h) (z)=0 with H^(z) = nz 2 j Q dk% A Hf(z, k*). 



(43) 



We assume that the sum rule holds in a strong sense, i.e., for pions and kaons separately and at the integrand level, 
for each value of z and kr- For pions, it follows that 

Hu {1/2) ,=-i^ (1/2) +v (44) 

l(u— >1T ) 1{U— >7T+) V ' 

The other u, d, d, unfavored Collins functions arc related to the above by isospin and charge symmetries, Eq. ([9]). Our 
strong interpretation of the Schafer-Teryaev sum rule together with Eq. (fT0|) (with D\ replaced by H^ 1 ^ 2 ^) implies 

<i?-) = -<^ + )=0. (45) 
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FIG. 7: Azimuthal asymmetry Ax%{z\,z%) for the production of two pions as a function of Z2 and integrated in bins of Zi at 
Q 2 = 110.7 GeV 2 . Dashed lines are obtained assuming Eq. (f37)l , solid lines assuming Eq. (|38]i . Note that the last zi bin in our 
calculation is narrower than in the corresponding experimental measurement. 



For kaons, the same considerations lead to the following assumptions 

<i 2 i-)=-<^ +)5 (47) 

rr-L(l/2) _ zx-L(l/2) _ n / 4 o\ 

In Fig. [7] we show the values of the pion azimuthal asymmetry for four different ranges of z\, as a function of z%- The 
dashed curves and solid curves are obtained respectively under the assumptions in Eq. ([37| and Eq. (f38|) . respectively. 
The upper curves exceed the data for the higher values, which either reflects the need for corrections to the linear 
approximation in Eq. (|40p . or more likely that assuming no evolution for the Collins function may be too severe an 
approximation. 

We calculated the corresponding KK asymmetry, Fig. and obtained even larger values, suggesting that there 
will be more dramatic effects in this accessible channel. 



V. CONCLUSIONS 



In this paper, we performed a new calculation of the Collins fragmentation function for u — > 7r + , along the lines 
of Refs. [n], [l2j but with some important differences: (i) we assumed the mass of the spectator is different from the 
mass of the fragmenting quark, (ii) we introduced a form factor at the hadron-quark vertex, (iii) we fitted the values 
of the model parameters to reproduce the unpolarized fragmentation function D\ at a scale Q\ — 0.4 GeV 2 . We 
compared the results of our model calculation to the available parametrizations of the Collins function [|| @] extracted 
from e + e~ annihilation and SIDIS data and found a reasonable agreement. We stressed the importance of critically 
considering different assumptions on the evolution of the Collins function with the energy scale. 

For the first time we presented an estimate of the Collins function for u — * K + and s — » K + . In particular, we 
found that the ratio H^ 1 ^ 2 * 1 / D\ for u — > K + is almost identical to that for u — > 7r + , while the ratio for s — > K + is 
about twice as big. 

Using the results of our model, we presented estimates for pion and kaon Collins asymmetries in e + e" annihilation at 
the BELLE experiment. In order to calculate the unfavored Collins functions we adopted the "strong interpretation" 
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FIG. 8: Azimuthal asymmetry A\2{z\, Z2) for the production of two kaons as a function of Z2 and integrated in bins of z\ at 
Q 2 = 110.7 GeV 2 . Dashed lines are obtained assuming Eq. (|37[) . solid lines assuming Eq. (|38[) . 



of the Schafer-Teryaev sum-rule 35]. Our results are in qualitative agreement with the available BELLE data on the 
pions, but large uncertainties arise from making different assumptions on the evolution of the Collins function as well 
as from determining the unfavored Collins fragmentation function. For the kaons, we predict the asymmetries to be 
larger than the pions. 
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